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We develop a theory for polymer translocation driven by a time-dependent force through an oscil¬ 
lating nanopore. To this end, we extend the iso-flux tension propagation theory (IFTP) [Sarabadani 
et al, J. Chem. Phys., 2014, 141, 214907] for such a setup. We assume that the external driving 
force in the pore has a component oscillating in time, and the flickering pore is similarly described 
by an oscillating term in the pore friction. In addition to numerically solving the model, we derive 
analytical approximations that are in good agreement with the numerical simulations. Our results 
show that by controlling either the force or pore oscillations, the translocation process can be either 
sped up or slowed down depending on the frequency of the oscillations and the characteristic time 
scale of the process. We also show that while in the low and high frequency limits the translocation 
time r follows the established scaling relation with respect to chain length No, in the intermediate 
frequency regime small periodic fluctuations can have drastic effects on the dynamical scaling. The 
results can be easily generalized for non-periodic oscillations and elucidate the role of time dependent 
forces and pore oscillations in driven polymer translocation. 


I. INTRODUCTION 

Polymer translocation through a nanopore, which has 
rapidly emerged as one of the most active research areas 
in biological and soft matter physics, plays an important 
role in many biological processes such as RNA transport 
through a nuclear membrane pore and a-hemolysin chan¬ 
nels [1], [2|, virus injection into cells and transportation of 
protein through membrane channels Q. It also has many 
technological applications such as rapid DNA sequencing 
Iffll , drug delivery 0j and gene therapy, and has been 
motivation for many theoretical and experimental stud¬ 
ies 0, @-[13 ■ So far most of the theoretical studies have 
been focused on the case of pore-driven translocation of 
a polymer chain through a nanopore with a constant ra¬ 
dius and by a constant driving force, which is inherently a 
far-from-equilibrium phenomenon [28l l30M38l ]. However, 
some theoretical and experimental work has also been 
done where the width of the pore changes during the 
translocation process. For example, the active translo¬ 
cation of a polymer chain through a flickering pore has 
been investigated by means of molecular dynamics simu¬ 
lations to show that translocation through the pore with 
an alternating width and sticky walls exhibits more effi¬ 
cient translocation as compared to the static pore (lit ]. 
There are some biological examples of this such as the 
nuclear pore complex, which plays an important role in 
nucleocytoplasmic transport in eukaryotes 03 ■ Another 
example is the exchange of molecules between mitochon¬ 
dria and the rest of cell which is controlled by the twin- 
pore protein translocase (TIM22 complex) in the inner 
membrane of the mitocondria 03 ] . On the experimental 
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side, it has been shown that the translocation of DNA 
through a nanochannel can be modulated by dynamically 
changing the cross section of an elastomeric nanochannel 
device by applying mechanical stress 02l - [44( . Moreover, 
the width of the nanopore can be tuned by using ther¬ 
mally driven nanoactuation of the polyNIPAM brushes 
inside the nanopore (45]. 

The external driving force can also be time-dependent 
during the translocation process [3, EfiUfiol ] . As an exam¬ 
ple, polymer translocation has been studied for a time- 
dependent alternating driving force by means of Langevin 
dynamics simulations. It has been shown that when the 
polymer-pore interaction is attractive, resonant activa¬ 
tion occurs at an optimal frequency of the alternating 
driving force 03 ■ There are some biological applications 
for the alternating driving force case such as translo¬ 
cation of linear and a-helical peptides throu gh a n a- 
hemolysin pore in the presence of an AC field 03, and 
monitoring the escape of the DNA from an a-hemolysin 
pore by using an alternating current signal 51]. More¬ 
over, using an alternating electric field in the nanopore 
has been suggested as a method for DNA sequencing [3]. 

On the theoretical side, over the last few years a con¬ 
sistent theory of driven translocation has been devel¬ 
oped based on the idea of tension propagation (TP) 0. 
Within the TP theory, the dynamics of driven transloca¬ 
tion can be described by considering a tension force that 
propagates along the backbone of the chain. The dynam¬ 
ics is dominated by the drag of the cis side of the poly¬ 
mer chain, leading to asymptotic scaling of the average 
translocation time t oc 7Vq + 1/ as a function of the chain 
length (number of monomers) N 0 03,01, S, Ell , where 
the Flory exponent v is exactly 3/4 in two and approxi¬ 
mately 0.59 in three dimensions, respectively 03| . In ad¬ 
dition, the theory reveals that the pore friction ?y p plays 
an important role for chains of finite length, leading to a 
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FIG. 1: (a) A schematic of the translocation process in the propagation stage, i.e. t < tt P ,T, for the trumpet regime. The 

time dependent external driving force f(t) acts only on the polymer bead in the pore towards the trans side. No is the length 
of polymer and translocation coordinate, s, is the number of beads that have already been translocated into the trans side. 
I + s is the number of beads influenced by the tension force on the cis side, which is less than No during the propagation stage. 
The location of the tension front is determined by R. d(t) is the pore width that is varying in time during the translocation 
process when the pore is flickering, (b) The translocation process for the trumpet regime during the post propagation stage 
that is the regime starting after which the tension front reaches the chain end. Panels (c) and (d) are the same as panel (a) 
but for the stem-flower and strong-stretching regimes, respectively. The propagation times for the trumpet, stem-flower and 
strong-stretching regimes are defined by f tp ,T, ft p ,sF and f tp ,ss, respectively. 


large correction-to-scaling term such that the transloca¬ 
tion time actually scales as t(Nq) = C\N P +V + C 2 ri p N 0 , 
where Ci and C 2 are constants [28l . [52j] . This scaling form 
has most recently been theoretically derived using the 
iso-flux Brownian dynamics tension pro pag ation (IFTP) 
theory of driven polymer translocation |25l |52| . 


In this paper, our main goal is to investigate the influ¬ 
ence of an alternating external driving force and a flick¬ 
ering pore on the dynamics of the driven translocation 
process. To this end we employ the IFTP theory to the¬ 
oretically study polymer translocation in the presence 
of a time dependent driving force through a flickering 
nanopore. Here, the time dependent driving force, f(t), 
is incorporated into the theory directly, while the effect of 
the flickering pore is implemented by a time dependent 
pore friction coefficient r] p (t). Using the IFTP theory 
modified in this manner allows us to explicitly derive the 
TP equations for the present setup. The formalism also 
allows us to investigate both the influence of thermal fluc¬ 
tuations from the solvent by introducing randomness into 
the effective driving force [27], [H, [52 and the influence of 
the distribution of the initial configurations of the chain 
on the cis side [H|. 


The outline of the paper is as follows: In Sec[TT]we show 
how to model pore-driven polymer translocation in the 
context of the iso-flux Brownian dynamics tension prop¬ 
agation formalism when the pore friction is time depen¬ 
dent and the driving force is alternating. Secs. IIII Al and 
IIIIB I present the finite-size scaling form for the transloca¬ 
tion time in the presence of an alternating driving force 
and time dependent pore friction, respectively. Secs. lIVI 
and m are devoted to the results on the average translo¬ 
cation time, the waiting time distribution, and time evo¬ 
lution of the translocation coordinate, respectively. Con¬ 
clusions and discussion are in Sec. lVIl 


II. MODEL 


For the sake of brevity, we use dimensionless units de¬ 
noted by tilde as Z = ZjZ u , with the units of length s u = 
a, time t u = rja 2 /(ksT), velocity v u = a/t u = fceT/^a), 
force f u = ksT/a 1 monomer flux <f> u = ksT / (rja 2 ), and 
friction r M = 77 , where a is the segment length, 77 is the 
solvent friction per monomer, ks is the Boltzmann con¬ 
stant, and T is the temperature of the system. 

To elucidate the parameter values chosen in this paper, 
we compare them with those in real physical systems and 
also with those used in molecular dynamics (AID) simula¬ 
tions in reduced Lennard-Jones (LJ) units [46l. l52i|. Here, 
parameters without tilde such as external driving force, 
F , pore friction, rj pi and time, t , are in LJ units. In MD 
simulations the scale of time, mass, and energy can be 
fixed by the LJ parameters m, a and e which are the mass 
and diameter of each bead, and the strength of the inter¬ 
action, respectively. We have chosen fcgT = 1.2 which is 
dimensionless and the time is scaled by = (ma 2 /e) 1 / 2 . 
In our model the interaction strength is 3.39 x 1(D 21 J 
at room temperature (T = 295 K), the mass of each 
bead is about 936 amu, and the size of each bead corre¬ 
sponds to the single-stranded DNA Kuhn length, which 
is a « 1.5 nm. Therefore, the LJ time scale is 32.1 ps. 
Here the thickness of the pore is set to a. By assuming 
the effective charge of 0.094e for each unit charge, three 
unit charges per bead [13, HH , and with a force scale of 
2.3 pN, an external driving force of F = 5 corresponds 
to a voltage of 375 mV across the pore fil ]. 

We use Brownian dynamics (BD) in the overdamped 
limit as the basic theoretical framework [23,111. Th e 
time evolution of the translocation coordinate s which is 
the chain length on the trans side is governed by the BD 
equation that can be written as 
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where T(t) is the effective friction, the first term in the 
right hand side is the entropic force and 7 ' is the surface 
exponent ( 7 ' ~ 0.95,« 0.69 for self-avoiding chains in 
two and three dimensions, respectively, and 7 ' = 0.5 for 
ideal chains), Nq is the chain length or the total number 
of beads in the chain (here the contour length of the 
polymer chain is L = aNo ), fit) is the external driving 
force, /tot is the total force, and the Gaussian white noise, 
C(Z), satisfies (£(£)) = 0 and (£(£)£(£')) = 2T(t)kBTS(t — 
t'). 

As both the cis side subchain and the pore friction 
contribute to the effective friction T(£), it can be writ¬ 
ten as r(f) = fjcisit) + Vp (£)■ For a fully flexible self¬ 
avoiding polymer chain, the dynamical contribution of 
the trans side to the dynamics has been shown to be in¬ 
significant [13, [H, [H, M E 3 and it can be absorbed into 
the pore friction fj p (t). Note that here the pore friction 
fj p (t) is assumed to be a function of time. We solve the 
dynamics of the cis side by using the TP equations. Fol¬ 
lowing the arguments of Rowghanian et al. [25[ we are 
able to derive the TP equations analytically. We assume 
that the monomer flux, <f> = ds/dt, is constant in space 
on the mobile domain of the cis side and through the 
pore, but evolves in time. The tension front which is the 
boundary between the immobile and mobile domains, is 
located at the distance x = —R(t) from the pore. The 
external driving force inside the mobile domain is me¬ 
diated by the backbone of the chain from the pore lo¬ 
cated at x = 0 to the last mobile monomer N which is 
located at the tension front. Using the force-balance re¬ 
lation for the differential element dx located between x 
and x + dx, the magnitude of the tension force at the 
distance x from the pore can be obtained. One can in¬ 
tegrate this force-balance relation over the distance from 
the entrance of the pore to x, to obtain the tension force 
at x as f(x,t) = f 0 -<t>(t)x, where f 0 = ftot~fjp(t) 4>(t) is 
the force at the entrance of the pore (see Appendix A of 
Ref. [H! for technical details). As the drag of all the pre¬ 
ceding monomers diminishes, the mediated force closer 
to the the tension front becomes smaller and it vanishes 
at the tension front. 

According to the blob theory, the shape of the mo¬ 
bile part of the chain can be classified into three differ¬ 
ent regimes. For a moderate external driving force, i.e. 
Nq” <C fo *C 1 , the monomer number density close to 
the pore is greater than unity, and the chain shape re¬ 
sembles a trumpet-like shape with the narrow end closer 
to the pore. This regime is called trumpet (TR) which is 
depicted in Figs. [TJi and b. For a stronger external driv¬ 
ing force, 1 <C /o -C Nq, the force can only straighten a 
small part of the chain which is called the stem, and the 
part following it is called the flower. This regime is called 
stem-flower (SF) (see Fig. [[]:). Finally, for a very strong 
driving force, i.e. fo Nq , the force is strong enough to 
straighten the whole moving part of the chain, which is 
called the strong-stretching (SS) regime (see Fig. [lji). In 
the TR and SF regimes, the location of the tension front 
is at the farthest blob from the pore, while in the latter 


one, i.e. SS, the tension front is located at the farthest 
mobile bead of the chain from the pore as can be seen in 

Fig. CD 

The monomer flux can be written as a function of the 
total force and linear size of the mobile part by inte¬ 
grating the force balance equation over the whole mobile 
part: 


4>(t) 


if tot (t) 

fj p {i) + R{i) 


( 2 ) 


The effective friction can be obtained as 


f(t) = R(t) + rjp(f), (3) 


by combining the definition of the flux, <j> = ds/dt, with 
Eqs. (HJ) and ©. 

Equations ©, m and (O together with the knowl¬ 
edge of R(t) determine the full solution for s. Therefore, 
the only equation remaining is the equation of motion for 
R(t) that must be derived separately for the propagation 
and post propagation stages. In the propagation stage, as 
shown in Fig. [TJi, c and d, the tension has not reached the 
last monomer of the chain, i.e. N = l + s < Nq, where 
N is the number of total beads that have been already 
influenced by the tension force and l is the number of 
dragged monomers in the cis side. Using the scaling rela¬ 
tion of the end-to-end distance of the self-avoiding chain, 
R = A V N”, where A v is a constant prefactor and N is 
the last monomer inside the mobile domain at the ten¬ 
sion front, together with the blob theory and N = l + s 
(see Appendix B of Ref. [52] for details) the equation of 
motion for the tension front in the propagation stage can 
be derived as 


R(t) 


vA„ R{i) * (jj(t) 
1 — vAh R(t)~~ 


(4) 


for the SS regime, and 


R(t) 


vAh RjiY^ 1 {(C a ,+G a .)[ftot(i)-4>(t)fjp(i)]+4>{i)} 

1 + is Ah R(t)~x~ £ a x </(£) 

(5) 


for the TR and SF regimes. Here the subscript ”a” in 
C a and g a stands for the trumpet regime as TR + and 
TR_ correspond to positive and negative values of <f>, 
respectively. For the stem-flower regime ”a” denotes SF. 

The time derivative of ftot is /tot- U a and Q a are given 
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by 


£tr+ 


£tr_ 


1?TR + 

£tr_ 

^SF 


!?SF 


fj p {t) + R(t) 

v P (t) 


[4>(t)R{t)] " 


4>(t) 




(2u - 1 )m 2 

: (j>(t) > 0, (6a) 


[ - <t>(t)R(t)\ 
[4>(t)R(t)] ~ 


fjp{t) + R(t) l (2v - 1 )4>{t) 2 
Vp {*) 

l 

4>{t) 
l 

4>(t) 


[-4>{t)R(t)] v 
cj){t) < 0 , (6b) 


: <j>(i) > 0, 


[ - cj>{t)R{t)] u : 4>{i) < 0, 


v — 1 


(2v - 1) [fjp(t) + 77(f)] </>(t) 2 

1 

4>{t) 


(6c) 

(6d) 

(be) 

(6f) 


In the post propagation stage which is schematically 
shown in Fig. ©) for the TR regime, all monomers on 
the cis side have already been affected by the tension 
force. Therefore, the correct closure relation for the post 
propagation stage is l+s = Nq. Therefore, one can derive 
the equation of motion for the tension front as 


R(t) = -m, ( 7 ) 


for the SS regime, while for the TR and SF regimes it 
can be written as 


R(i) 


(^a + Ga) [ftotjt) - 4 >(j)fj p (i)] + 4 >{t) 
4>(t) X £ a 


( 8 ) 


For more details, the interested reader is referred to the 
Appendix B of Ref.l52l. 

Eqs. (Ill) , {1771) , m , (141) and © must be solved self- 
consistently to find the solution for the model in the 
propagation stage, while in the post propagation stage, 
one must solve the set of Eqs. ©, ©, ©, © and ©. 


III. SCALING OF TRANSLOCATION TIME 
FOR SS REGIME 

A. Static pore and oscillating force 

To obtain analytical results, we assume that only the 
external driving force has contribution to the total force 
in the BD equation ©, i.e. we ignore both the entropic 
and thermal noise terms here. We choose the external 
driving force as a combination of a constant force F and 
an oscillatory term Af (t) as 

/tot = f(t) =F + A f (t), (9) 



FIG. 2: Analytic approximation for the normalized translo¬ 
cation time, r/f s tat, as a function of normalized force fre¬ 
quency, <J//u>stat for various values of initial phases of ■)/>/ = 
0, 7t/2,7t and 37r/2 in the SS regime. Here, the pore is static 
with constant pore friction r/po = 3.5, the chain length is 
No ~ 128, and the time dependent external driving force is 
f(t) = F + 5/ sin(oj ft + ipf) , where F = 25 and a/ = 2.5. The 
green diamonds show the normalized translocation time for 
the full IFTP model which includes the entropic force term. 
The magenta dash-dash-dotted curve presents the normalized 
translocation time in the small amplitude approximation. 


where Af(i) = df sin(w/f + ipf), and df, Qf and ifjf are 
the amplitude, frequency and initial phase of the force, 
respectively. For simplicity, we first consider a static pore 
with a constant pore friction, i.e. = fj p o. Then 

Eq. © reduces to (f>(t) = [F + Af(t)]/[R(t) + 77 p0 ]. In 
the strong-stretching regime considered here, the number 
of mobile monomers on the cis side is given by l$s = R- 
Knowing the external driving force as a function of time 
and pore friction together with the conservation of mass, 
N = s + l, the propagation time r t p can be solved by 
integration of N from 0 to Nq. The result is 




R(N)dN + fj p0 N 0 


Af tp , 


( 10 ) 


where 


A r tp =j 


^R 2 (N 0 )+ijp 0 R{N 0 )+ j tP A f (t)dt 


. ( 11 ) 


In the post propagation stage, one sets the condition 
dN/dt = 0 and integrates R from R{Nq) to 0 to obtain 
the post-propagation time, f pp , as 


F 


-i? 2 (./Vo)-l-?) p o-R(-/Vo)— / Af(t)dt 


’ r tp 


( 12 ) 
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Finally, the total translocation time is then given by 


t = f t p + f p 


1 

F 


r N o 


R(N)dN + rj p0 N 0 — / A f (t)dt 

Jo 


— Tstat ~ / Af(t)dt. 
r Jo 


(13) 


where 


f stat = . N , 0 1+ ^ + Mn 0 , 

(1 + v)F 0 F 


(14) 


is the translocation time in the absence of the oscillatory 
term in the external driving force. Here we have used the 
same value of the amplitude A„ = 1.15 as in Ref.fiisl. 

In Fig. [2] we plot the translocation time from Eq. (fl3l) 
normalized with Eq. m, as a function of the normal¬ 
ized frequency uif/Q B tat for various values of ipf in the 
SS regime (F = 25 and a/ = 2.5). For comparison, we 
also show results from the full deterministic IFTP model 
(green diamonds) where the entropic force term has been 
included. It can be seen that the agreement with the an¬ 
alytic approximation in the SS regime is very good for a 
wide range of frequencies. 

Using the explicit form of Af(t) = a/sin(w/f + i/jf) 
we can now examine Eq. m in different limits. First, 
in the high frequency limit, i.e. ujf t — , the inte¬ 
gral fyAffydt simply vanishes and the total transloca¬ 
tion time in the presence of the oscillatory term is the 
same as the translocation time in the absence of it, i.e. 
f = fstat, as can also be seen in the high-frequency limit 
of Fig. [2] 

In the opposite limit of a very small frequency, u>f <C 
t , the integral gives 


lj%m = -i 


a f 


FCjj 


df sin {id ft + ipf)dt 
cos ('*/’/) — COS (Cjf T +1pf ) 


.(15) 


The expansion of cos(w/f + ipf) to first order yields 
cos(ojfT+il>f) ~ cos(ijjf) — diff and thus 



FIG. 3: The normalized translocation time, f/f, stat, as a 
function of normalized force frequency, ujf/u) s tat for various 
values of initial phases of ijjf = 0, n/2, n and 3ty/2. Here, the 
pore is static with constant pore friction r/ p o = 3.5, the chain 
length is No = 128, and the time dependent external driv¬ 
ing force is f(t) = F + df sin(a jft + tpf), where F = 5.0 and 
af = 0.5. The orange circles show the normalized transloca¬ 
tion time for the full model where includes thermal noise and 
initial shape distribution of the chain. 


replaced by f s t a t- The translocation time can then be 
written as 


7" — 't’stat 


1 + 


o/ C0S ( 27r 3^T + V>/) - cos(V>/) 


2t tF 




• (17) 


This equation predicts an oscillatory dependence of the 
translocation time on the driving frequency. In Fig. [5] we 
plot the small amplitude expansion result for the case 
where the amplitude ratio af/F = 0.1 (magenta dash- 
dash-dotted line). For this relatively small ratio, the 
agreement with the full solution of the model (orange 
curve) is good for a wide range of frequencies. 


B. Flickering pore and oscillating force 


[l + d f sin (ip f )/F] 

As can be seen from Eq. (1131) , the scaling form for the 
average translocation time is the same for both static 
and high frequency limits. In the low frequency limit, 
the scaling with respect to the chain length is also the 
same, but the scaling form with respect to the external 
driving force becomes f ~ [F + df sin(^/)] _1 . 

Finally, if the amplitude of the oscillating term, df, is 
much smaller than the static force, F, then we can use 
a small amplitude approximation, where the value of the 
upper limit in the integral over time in Eq. (1131) can be 


We consider next the most general case of an oscil¬ 
lating force and a flickering pore in the SS regime. We 
again neglect both the entropic force and thermal noise, 
and take the oscillatory force to be the same as in the 
previous section. The flickering pore has the following 
time dependent pore friction: 

rj p {t) = f] p o + A p (t), (18) 

where A p (t) = a p sin(u;pt + ^> p ) and a p , d> p and ip P are 
the amplitude, frequency and initial phase of the pore 
friction, respectively. Then, Eq. © reduces to <j>{t) = 
[F + Af(t)]/[R(t) + 7? p0 + A p (t)]. As in Sec. lIII A1 we 
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FIG. 4: The normalized translocation time, r/fstat, as a func¬ 
tion of normalized frequency of the flickering pore, w p /w 3 tat, 
for various values of initial phases of ip p = 0,ir/2,n and 
37 t/ 2. Here, the driving force is constant, F = 5.0, the chain 
length is No = 128, and the time dependent pore friction is 
rj P {t) = Vpo+o-p sin(tUpf-l-V’p), where r / p o = 3.5 and a p = 1.275. 
The arrows indicate which way the pore moves at the begin¬ 
ning of the drive, due to the initial phase ip p . 


solve the propagation time by integration of N from 0 to 
Nq, and the post-propagation time by integration of R 
from R(Nq) to 0. The total translocation time is then 
given by 

1 _ _ 1 r N °~ _ 

r = T sta t--H / A f (t)dt+ / A p (t)dN. (19) 
r Jo T Jo 


Using the explicit form of A p (t) = a p sin(u} p t + ip p ), at 
the high frequency limit, w p t _1 , and for very long 
chains the integral J r ^°A p (t)dN vanishes and the whole 
translocation time in the presence of the oscillatory force 
is the same as the translocation time with a constant pore 
friction fj p o- We note, however, that for short chains the 
integral f^°A p (t)dN can be non-zero due to the interac¬ 
tion of the polymer beads with the pore [H]. 

In the opposite limit of a low frequency the integral 
gives 


1 rN o 


r N o 


F 


Therefore, 


A p {t)dN = — / a p sin(tD p t + ip p )dN 


F 

g p sin('0 p ) 
F 


N 0 . 


F 


r No 


R(N)dN + {? 7 P o + a p sm(ip p )}No 


( 20 ) 


— / Af{t)dt 


— Jgtat T 


d p sm(i/j p )No - Af(t)dt . 
Jo 

( 21 ) 


As we can see in Eq. m, at the low frequency limit 
the pore behaves like a static pore with pore friction of 
Vpo + d p sin(^p). 

Equations m and m show that the scaling form 
for the average translocation time for a very long chain 
at the high and also at the low pore frequency limits is 
the same as in the case where only the driving force is 
alternating and the pore is static (as discussed in the 
previous subsection). For a short chain at both the high 
and the low pore frequency limits the scaling may be 
affected by the value of the integral (1/E) J^°A p (t)dN. 

When the driving force is constant, if the amplitude of 
the oscillating term, d p , is much smaller than the static 
pore friction, f] p o, then we can also use here a small am¬ 
plitude approximation, where the value of the time in the 
second integral in Eq. (fliil) must be replaced by t(N(fj p0 )). 
Then, for A p (i) = d p sin(d} p t + ip p ), Eq. (TTUl) can be cast 
into 


r R(N 0 ) 


— fstat T 


uFAl /u Jo 


dR R 


sinjwp t(R) + i> p ], 

( 22 ) 


where t(R) = 
result of Eq. 


- I±£ 

R - 


(I+ „> - * + (5£ - %p° 


The 


has been checked and it matches very 
well with the result of the full deterministic model in a 
wide range of frequencies. 


IV. AVERAGE TRANSLOCATION TIME 

Next we consider the average translocation time, f, in 
all the different regimes using the full IFTP theory. To 
this end, we numerically solve the IFTP Eqs. ©, ©, 
©, ©, ©, 0 and © to obtain the total translocation 
time. We study the behavior of the translocation time 
first for a time-dependent force, time-dependent pore fric¬ 
tion, and finally when both the force and the pore friction 
are functions of time. 

In Fig. [3] the normalized translocation time, f/f s t a t, 
has been plotted as a function of normalized frequency, 
w//w s tat for various values of the initial phases of ipf = 
0 ,7 t/2 ,7 t and 3n/2. The pore is static with the pore 
friction r ] p o = 3.5, chain length is Nq = 128, and 
the time dependent external driving force is given by 
/(/) = F + df sin(Cjft + ipf), where F = 5.0 and a/ = 0.5. 
Here, tDstat = 27r/f s t a t is the frequency related to the 
whole translocation process when the pore is not flick¬ 
ering, i.e. fj p (t) = fj p o and the external driving force 
is constant, i.e. f{t) = F. For ipf = 0 (black solid 
curve) at low frequencies the magnitude of the translo¬ 
cation time is the same as of the static force case. At 
this low frequency limit the translocation process occurs 
during the first half of the period of the oscillating force. 
During this first half of the cycle the value of the force in¬ 
creases from its static value to its maximum and back to 
the static value at 0Jf/u> s t a t = 0.5, being always greater 
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FIG. 5: (a) The normalized translocation time, f/rstat, as a function of either normalized pore frequency, Wp/aistat, or force 
frequency, w//w a tat, for mixed initial phases tjjf = 7r/2 and ip p = On/2, ipf = On/2 and ip p = 7r/2, ipf = 7r/2 and ip p = n, 
and ifjf — 7t and ip p = 7r/2. Here, the pore is flickering and the pore friction is given by f/ P (t) = fj p o + a p sin(w p t + ip P ), where 
r/ p0 = 3.5 and a p = 1.275 while the oscillating external driving force is f(t ) = F + df sin (uift + ipf), where F = 5.0 and df = 0.5. 
(b) The normalized translocation time as a function of either c5 p /aj a tat or u)f/Q a ta.t, but for the identical initial pore and force 
phases ipf = i/j p = 0,n/2,n and On/2, with the same parameters as of (a). The arrows show the direction of the motion of a 
point on the pore wall at the starting of the translocation process. For both panels (a) and (b) the chain length is No = 128 
and 0 J p / u)stat bjf /c^stat - 


than its static value. Therefore, the translocation time 
gradually decreases for small frequencies, and has a min¬ 
imum at Cjf/Q a tat = 0.5. For frequencies higher than 
Cjf/Cj stat = 0.5, the chain starts experiencing the second 
half of the cycle where the value of the force is always 
smaller than its static limit, and thus there is a first 
maximum at a}//(D s tat = 1- For higher frequencies, i.e. 
uj f /cDstat > 1, the chain again experiences the following 
half of the cycle, i.e. Tf < t < 3X//2 (Tf = 2ir/ujf), 
where again the value of the force is greater than its static 
limit and the translocation time increases. With increas¬ 
ing frequency the translocation time displays oscillations 
between minima and maxima with a decreasing ampli¬ 
tude upon approaching the high-frequency limit, where 
the rapidly oscillating force component averages to zero, 
as discussed in Sec. nm For comparison, we have also 
solved the full stochastic IFTP model for the initial phase 
ipf = 0, with results shown in orange circles in Fig. [3] In 
this case we take into account both thermal noise and 
the initial shape distribution of the chain as explained in 
detail in Ref.[52. The results are in very good agreement 
with the deterministically solved equations. 

Next, we consider the case of a flickering pore with 
a constant drive. In Fig. [4] the normalized transloca¬ 
tion time, f/f s tat, has been plotted as a function of 
normalized frequency of the flickering pore, w p /w s tat, 
for a constant driving force, F = 5.0, chain length of 
No = 128, and time dependent pore friction fj p (t) = 
rj p0 + d p sin(w p t + i/j p ) with static pore friction rj p o = 3.5 
and amplitude a p = 1.275, for various values of the 
initial phases of ip p = 0,7r/2,7r and 3n/2. As we can 


see for ipp = 0, in the low frequency limit the value of 
the translocation time reduces to its static value, where 
t/ T s tat = 1. As the pore oscillation frequency increases 
the value of the pore friction also increases, which means 
that the pore width, d(t), decreases (see Fig.[T|a,). There¬ 
fore, the effective interaction between the pore and the 
beads becomes stronger and the result is that the value of 
the translocation time grows and reaches its maximum at 
w p /w s tat = 0.5. At this frequency the translocation pro¬ 
cess occurs during the first half period of the pore motion. 
During this first half of the cycle the radius of the pore 
decreases from its static value to its minimum and back 
to the static limit at w p /w s t a t = 0.5, being always smaller 
than the radius of the static pore. Therefore, the translo¬ 
cation time gradually increases for low frequencies, and 
then it has a maximum at c5 p /u; s tat = 0.5. For frequen¬ 
cies higher than w p /w s t a t = 0.5, the chain experiences the 
second half of the cycle where the radius of the pore is 
larger than the radius of the static pore, and thus there 
is a strong minimum at w p /w s tat = 1- For even higher 
frequencies, i.e. <I> P /(D S tat > 1, the chain experiences the 
next half of the cycle, i.e. T p < t < 3T p /2 (T p = 2 tt/u p ), 
where the radius of the pore is again smaller than its 
static value and the translocation time increases. This 
alternation of high and low pore friction creates succes¬ 
sive minima and maxima in a way similar to the case 
of the time-dependent drive, except that increasing the 
pore friction leads to increasing translocation time, while 
increasing the drive leads to more rapid translocation. 

Finally, we consider the combination of an oscillating 
drive and pore. The fact that increasing the drive or 













FIG. 6: The effective translocation exponent, a(No), as a 
function of the chain length, No, for two different values of ini¬ 
tial phases ipf = 0 (red dashed curve) and n/2 (green dashed- 
dotted curve). The black solid circle-line curve shows the 
static case where both the pore friction and the external driv¬ 
ing force are constant with the values r/ p o = 3.5 and F = 5.0. 
Red and green curves are for the cases where the pore fric¬ 
tion is constant, r/ p o = 3.5, but the external driving force is 
time-dependent and given by /(f) = F + 2/sin(a)/f + ipf), 
with F = 5.0, a/ = 0.5 and fixed uif = 0.5o> s tat (No = 1000) = 
9.5 x 10~ 4 . The horizontal blue solid line is the asymptotic 
value of the translocation exponent which is 1 + v in the limit 
of very long chains. 


increasing friction work in the opposite directions with 
respect to the translocation time suggests that it could 
be minimized by out-of-phase combination of the two 
forces. In Fig.[5jr we present the normalized transloca¬ 
tion time as a function of either the normalized pore 
frequency or force frequency, for mixed initial phases 
ipf = 7t/ 2 and ip p = 37r/2, ipf = 37r/2 and ip p = n/2, 
ipf = 7t/2 and ip p = n, and ipf = n and ip p = n/2. 
Here, the pore is flickering and the pore friction is given 
by Vp(t) = Vpo + «p sin(d}pt + ipp), where r/po = 3.5 and 
a p = 1.275 while the oscillating external driving force is 
/(f) = F + af sin(d )ft + ip /) where F = 5.0 and a/ = 0.5. 
The pore frequency is equal to the force frequency dur¬ 
ing the translocation process here. For fixed values of 
F, a,f, rjpQ and a p , the minimum value of the normal¬ 
ized translocation time, T m i n /f s t a t, occurs at very low 
frequency limit when the force phase is ipf = n/2 and 
pore phase is ip p = 3n/2, because the driving force has 
its maximum value and the pore friction has its min¬ 
imum value. Then the normalized translocation time 
gets its maximum value, f max /f s t a t, when ipf = 3n/2 
and pore phase is ip p = n/2, because the driving force 
has its minimum value while the pore friction has its 
maximum value. All other values of the normalized 
translocation times, r/f s t a t, for the other values of the 
force and the pore phases at different force or pore fre¬ 
quencies are between these two extremum values, i.e. 


Arliri /"Gtat A t/t s t a t 7"M a x/^”st a t ■ 

In Fig. [5] b we show the same quantity as a function 
of either (I> p /w s t a t or Cbf/uj s t a t, but the initial pore and 
force phases are now identical, and />/ = ip v = Q,n/2,n 
and 3n/2, with the same parameters as of Fig.[5]i. Both 
the pore friction and the driving force are time depen¬ 
dent with the same forms and the same parameters as of 
Fig. [5ja. As expected, the two effects now work against 
each other, and there’s an almost complete cancellation 
in the low and high frequency limits. As the effect of time 
dependent driving force and pore friction are dynamically 
coupled to each other, it is very difficult to find parame¬ 
ters in such a way that at low frequency limit they exactly 
cancel the effect of each other. For intermediate frequen¬ 
cies there are some oscillations due to the fact that the 
total effective friction, Eq. ([3]), is determined dynamically 
by tension propagation on the cis side chain, and cannot 
be completely controlled externally. In the end, adding a 
time-dependent component to the pore friction and the 
drive is fundamentally different. 

Finally, it is interesting to see how the (effec¬ 
tive) translocation exponent a, defined as a(No) = 
d Inf/din No, is affected by the presence of an alternat¬ 
ing external driving force. For a static driving force, the 
correction-to-scaling term due to the pore friction is large 
here and a slowly approaches the value 1+zz s=s 1.588 with 
increasing chain length [H US E3, as shown in Fig .[6] 
For the oscillating force, f(t) = F + d/sin(w/t + ipf) 
with F = 5.0, af = 0.5, constant pore friction ? 7 p o = 3.5, 
and fixed w/ = 0.5a} sta t(lVo = 1000), we see a very 
strong oscillatory behavior of a. The chosen frequency 
Cbf = 0.5£f s t a t(-/Vo = 1000) = 9.5 x 10 -4 is moderate com¬ 
pared to the translocation time of the chain of length 
A/) = 1000, low for a short chain such as Nq = 40, and 
very high for a long chain such as A^o = 10 5 . By fix¬ 
ing the frequency when the chain length is changed, the 
limits of low, moderate and high frequencies are all ex¬ 
plored. At the very short and very long chain limits the 
translocation exponent is the same as for constant driv¬ 
ing and independent of the value of the initial phase ipf. 
However, close to Nq = 1000, the translocation exponent 
oscillates due to the oscillation of the translocation time 
(see Fig.[3|). It is worth noting that while the amplitude of 
the exponent’s oscillation depends on the ratio of af/F, 
even with a moderate value of the ratio ofa//F = 0.1 the 
exponent varies between 1.40 and 1.74, going well beyond 
1 + v. In addition, the oscillation amplitude decreases 
towards long chains, as the system effectively enters the 
high frequency limit and the average driving force ap¬ 
proaches the static force value. 


V. DYNAMICS OF THE TRANSLOCATION 
PROCESS 

To understand the effect of the oscillating pore and 
force on the translocation process, it is necessary to ex¬ 
amine also the dynamics of the process and not just the 
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FIG. 7: (a) The waiting time, w(s), as a function of the translocation coordinate, s, for various values of initial force phases 
ipf = 0, 7r/2,7 r and 3n/2. Here, the alternating external driving force is given by f[t) = F + 5/ sin(oi/t + />/), where F = 5 and 
a/ = 0.5, the pore friction is constant, r/ p o = 3.5, and the force frequency is Cbf = 0. Panels (b), (c) and (d) are the same as 
(a) but for different values of normalized force frequencies uif /w s tatic = 0.2,1 and 4, respectively. In all panels the chain length 
is No = 128, and the black solid curve has been plotted for completely static case where the external driving force is constant, 
F = 5, and the constant pore friction is %,o = 3.5. 


resulting translocation time. Often the translocation pro¬ 
cess has been characterized by calculating the time evo¬ 
lution of the translocation coordinate s. However, a more 
sensitive measure is the monomer waiting time distribu¬ 
tion that measures the time a monomer spends inside the 
pore. The tension propagation (increasing waiting time) 
and post propagation (decreasing waiting time) stages 
are also easy to distinguish from the distribution. For 
the oscillating force and static pore, the waiting time dis¬ 
tributions are shown in Fig.0 For the oscillating pore, 
the waiting times are not shown, but they are analogous 
to the oscillating force and consistent with the results 
discussed in the previous section. 

In Fig.[71i the waiting time, w(s), is plotted as a func¬ 
tion of the translocation coordinate, s, for various val¬ 
ues of initial force phases ipf = 0,7r/2,7r and 3tt/2. 
Here, the alternating external driving force is given by 
/(f) = F + af sin(w/f + ijjf), where F = 5 and a/ = 0.5, 


the pore friction is constant, r) p o = 3.5, and the force 
frequency is w/ = 0. The chain length is No = 128, and 
the black solid curve has been plotted for the completely 
static case, where the external driving force is constant, 
F = 5, and the constant pore friction is r/ p o = 3.5. With 
the initial phases ipf = 0 and 7r, the waiting time dis¬ 
tribution curves collapse on the waiting time distribu¬ 
tion curve for the static case, and as a result the corre¬ 
sponding translocation times for these two initial phases 
are exactly the same as it can be seen in Fig. [3] With 
ijjf = 7t/2 the value of the force at w/ = 0 is greater 
than its value for the static case and therefore the wait¬ 
ing time is consistently below the static force curve and 
the overall translocation time is smaller. For ijjf = 37t/2 
the converse is true. 

In Figs. 03, c and d the waiting time distributions are 
shown for low, intermediate and high frequencies. In the 
low frequency regime, the waiting time gradually devi- 
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ates from the static force distribution, while for the high 
frequency the distribution oscillates between the large 
and small force values. As the translocation time is the 
integral of the waiting time, it is clear from the high 
frequency distribution (Fig.®) why the effect of the os¬ 
cillating force is averaged out. On the other hand, for 
the intermediate frequencies the waiting time and thus 
the translocation time can be higher or lower depending 
on the initial phase. 

In addition to the waiting time, we have also solved 
for the time evolution of the translocation coordinate s. 
In this case, we have used the full stochastic model in¬ 
cluding the fluctuations due to thermal noise and initial 
configuration such that we can extract the variance of s 
in addition to the mean value. In addition, we have also 
increased the driving force and the amplitude to make 
the oscillations more visible. The simulation parameters 
are thus fit) = F + a/sin(w/t + iff), where F = 7, 
a/ = 3 , iff = 7T, pore friction 7y p o = 3.5, and chain length 
No = 128. The results for the mean value (s(f)) and the 
variance (<5s 1 2 3 4 5 6 7 8 9 10 11 (i)) = (s 2 (t)) — (s(t)) 2 are shown in Figs® 
c. The oscillations in the mean value of s(t) are visible for 
the intermediate frequencies similar to the waiting time 
distribution. 


VI. CONCLUSIONS 

We have extended the recently developed iso-flux ten¬ 
sion propagation theory to polymer translocation driven 
by a time-dependent external force and through an os¬ 
cillating nanopore. Specifically, we take into account a 
periodically oscillating term in the external driving force 
and describe the flickering of the pore with a similar term 
in the pore friction. We have derived analytical approx¬ 
imations for the translocation time in the low and high 
frequency limits, and also in the intermediate frequency 


regime for small oscillation amplitudes. The analytical 
expressions can be used to interpret the numerical re¬ 
sults obtained from solving the full model. The theory 
can also be easily generalized for other types of oscilla¬ 
tions, including non-periodic ones. 

We have shown that in the low and high frequency lim¬ 
its, the dynamics of the translocation process reduces to 
the non-oscillating cases; the translocation time is found 
to scale similarly to the static force and pore simulations, 
following the finite size scaling relation f = c\Nq + 1 ' + 
c 2 r] p No. However, in the intermediate frequency regime, 
determined by the characteristic time scale of the process, 
the effective scaling exponent, a{No ) = dinf/din No, is 
shown to be extremely sensitive to small periodic fluc¬ 
tuations. The results indicate that such periodic fluc¬ 
tuations can indeed be used to either speed up or slow 
down the translocation process. Additionally, it is pro¬ 
posed that such control can be achieved either by tuning 
the oscillations in the driving force or in the pore char¬ 
acteristics. On the other hand, it is plausible that such 
oscillations may have significant influence on the experi¬ 
mentally measured translocation times and scaling expo¬ 
nents. 

Our results shed light on the role that small oscillations 
in the external driving force and in the pore characteris¬ 
tics have on the dynamics of polymer translocation. The 
developed theory and the results may help in understand¬ 
ing both biological functions and experiments, which are 
being done on increasingly complex nanopore systems. 


Acknowledgments 

This work was supported by the Academy of Finland 
through its Centres of Excellence Program (2012-2017) 
under Project No. 915804. 


[1] H. Salman, D. Zbaida, Y. Rabin, D. Chatenay and M. 
Elbaum, Proc. Nati. Acad. Sci. 98 7247 (2001). 

[2] J. J. Kasianowicz, E. Brandin, D. Branton and D. W. 
Deamer, Proc. Nati. Acad. Sci. 93 13770 (1996). 

[3] B. Alberts, A. Johnson, J. Lewis, M. Raff, K. Roberts 
and P. Walter, Molecular Biology of the Cell (Garland 
Science, 2002). 

[4] A. Meller, J. Phys.: Condens. Matter, 15 R581 (2003). 

[5] A. Meller, L. Nivon and D. Branton, Phys. Rev. Lett. 86 
3435 (2001). 

[6] K. Luo, T. Ala-Nissila, S.-C. Ying and A. Bhattacharya, 
Phys. Rev. Lett. 100 , 058101 (2008). 

[7] G. Sigalov, J. Comer, G. Timp and A. Aksimentiev, Nano 
Lett. 8, 56 (2008). 

[8] J. A. Cohen, A. Chaudhuri. and R. Golestanian, Phys. 
Rev. X 2, 021002 (2012). 

[9] W. Sung and P.J. Park, Phys. Rev. Lett. 77, 783 (1996). 

[10] M. Muthukumar, J. Chem. Phys. Ill, 10371 (1999). 

[11] J. Chuang, Y. Kantor and M. Kardar, Phys. Rev. E 65, 


011802 (2001). 

[12] R. Metzler and J. Klafter, Biophys. J. 85, 2776 (2003). 

[13] Y. Kantor and M. Kardar, Phys. Rev. E 69, 021806 
(2004). 

[14] A.J. Storm et al, Nano Lett. 5, 1193 (2005). 

[15] A. Yu. Grosberg, S. Nechaev, M. Tamm and O. Vasilyev, 
Phys. Rev. Lett. 96, 228105 (2006). 

[16] J.L.A. Dubbeldam, A. Milchev, V.G. Rostiashvili and 
T.A. Vilgis, Europhysics Lett. 79, 18002 (2007). 

[17] T. Sakaue, Phys. Rev. E 76, 021803 (2007). 

[18] T. Sakaue, in Proceedings of the 5th Workshop on Com¬ 
plex Systems , AIP CP 982, 508 (2008). 

[19] K. Luo, S.T.T. Ollila, I. Huopaniemi, T. Ala-Nissila, 
P. Pomorski, M. Karttunen, S.-C. Ying and A. Bhat¬ 
tacharya, Phys. Rev. E 78 050901 (R) (2008). 

[20] D. Branton, D.W. Deamer, A. Marziali et al, Nature 
Biotech. 26, 1146 (2008). 

[21] K. Luo, T. Ala-Nissila, S.-C. Ying and R. Metzler, Eu- 
rophys. Lett. 88, 68006 (2009). 



11 





FIG. 8: (a) The averaged translocation coordinate, s(t), as a function of time, t, for the static case with constant external 
driving force F = 7 and constant pore friction r/ p o = 3.5 (solid black curve), and also when the external driving force is 
alternating with /(f) = F + a/sin(d)/f + ipf) where F = 7, af = 3 and ipf = n, and the constant pore friction is r/ p o = 3.5, 
for various values of force frequencies 6)f /6) s ta.t = 0.5,1,2 and 4 (red dashed, green dashed-dotted, blue dashed-dotted-dotted 
and pink dashed-dashed-dotted curves, respectively). Panels (b) and (c) present the variance ( Ss 2 (t )) = ( s 2 (t )) — (s(t)) 2 as a 
function of time for the same parameters as of panel (a) in normal-normal and log-log scales, respectively. The chain length is 
No = 128 in all panels. 


[22] E. Schadt, S. Turner and A. Kasarskis, Hum. Mol. Gen. 
19, R227 (2010). 

[23] T. Sakaue, Phys. Rev. E 81, 041808 (2010). 

[24] A. Milchev, J. Phys.: Condens. Matter 23, 103101 

( 2011 ). 

[251 P. Rowghanian and A. Y. Grosberg, J. Phys. Chem. B 
115, 14127 (2011). 

[261 M. Muthukumar, Polymer Translocation (Taylor and 
Francis, 2011). 

[27] T. Ikonen, A. Bhattacharya, T. Ala-Nissila and W. Sung, 
J. Chem. Phys. 137, 085101 (2012). 

[28] T. Ikonen, A. Bhattacharya, T. Ala-Nissila, W. Sung, 
Europhys. Lett 103, 38001 (2013). 

[29] V. V. Palyulin, T. Ala-Nissila and R. Metzler, Soft Mat¬ 
ter, 10, 9016 (2014). 

[30] A. Bhattacharya, W.H. Morrison, K. Luo, T. Ala-Nissila, 

S. -C. Ying, A. Milchev and K. Binder, Eur. Phys. J. E 
29, 423 (2009). 

[31] V. Lehtola, R.P. Linna and K. Kaski, Europhys. Lett. 
85, 58006 (2009). 

[32] A. Bhattacharya and K. Binder, Phys. Rev. E 81, 041804 

( 2010 ). 

[33] V.V. Lehtola, K. Kaski, R.P. Linna, Phys. Rev. E 82, 
031908 (2010). 

[34] J. L. A. Dubbeldam, V. G. Rostiashvili, A. Milchev and 

T. A. Vilgis, Phys. Rev. E 85, 041801 (2012). 

[35] T. Saito and T. Sakaue, Eur. Phys. J. E 34, 135 (2012). 

[36] T. Saito and T. Sakaue, Phys. Rev. E 85, 061803 (2012). 

[37] T. Saito and T. Sakaue, arXiv: 1205.3861 (2012). 

[38] T. Ikonen, A. Bhattacharya, T. Ala-Nissila and W. Sung, 
Phys. Rev. E 85, 051803 (2012). 

[39] J. A. Cohen, A. Chaudhuri. and R. Golestanian, Phys. 
Rev. Lett. 107, 238102 (2011). 

[40] J. Yamada, J. L. Phillips, S. Patel et al, Mol. Cell. Pro- 
teomics 9, 2205 (2010). 


[41] P. Rehling, K. Brander and N. Pfanner, Nat. Rev. Mol. 
Cell Biol. 5, 519 (2004). 

[42] D. Huh, K. L. Mills, X. Zhu et al. , Nat. Mat. 6, 424 
(2007). 

[43] E. Angeli, C. Manneschi, L. Repetto et al, Lab. on Chip 
11, 2625 (2011). 

[44] P. Fanzio, C. Manneschi, E. Angeli et al., Sci. Rep. 2, 
791 (2012). 

[45] B. Yameen, M. Ali, R. Neumann et al, Small 5, 1287 
(2009). 

[46] T. Ikonen, J. Shin, W. Sung and T. Ala-Nissila, J. Chem. 
Phys. 136, 205104 (2012). 

[47] J. A. Cohen, A. Chaudhuri. and R. Golestanian, J. Chem. 
Phys. 137, 204911 (2012). 

[48] A. Fiasconaro, J. J. Mazo and F. Falo, Phys. Rev. E 91, 
022113 (2015). 

[49] R. I. Stefureac, A. Kachayev and J. S. Lee, Chem. Com- 
mun. 48, 1928 (2012). 

[50] M. Bates, M. Burns and A. Meller, Biophys. J. 84, 2366 
(2003). 

[51] D. K. Lathrop, E. N. Ervin, G. A. Barrall et al, J. Am. 
Chem. Soc. 132, 1878 (2010). 

[52] J. Sarabadani, T. Ikonen and T. Ala-Nissila, J. Chem. 
Phys. 141, 214907 (2014). 

[53] M. Doi and S. F. Edwards, The Theory of Polymer Dy¬ 
namics (Oxford University Press, 1986). 

[54] A. F. Sauer-Budge, J. A. Nyamwanda, D. K. Lubensky 
and D. Branton, Phys. Rev. Lett. 90, 238101 (2003). 

[55] J. Mathe, H. Visram, V. Viasnoff, Y. Rabin and A. 
Meller, Biophys. J. 87, 3205 (2004). 

[56] J. L. A. Dubbeldam, V. G. Rostiashvili and T. A. Vilgis, 
J. Chem. Phys. 141, 124112 (2014). 

[57] P. M. Suhonen, K. Kaski and R. Linna, arXiv: 1405.0920 
(2014). 











